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We extend the Berry-phase concept of polarization to insulators having a non-zero value of the 
Chern invariant. The generalization to such Chern insulators requires special care because of the 
partial occupation of chiral edge states. We show how the integrated bulk current arising from an 
adiabatic evolution can be related to a difference of bulk polarizations. We also show how the surface 
charge can be related to the bulk polarization, but only with a knowledge of the wavevector at which 
the occupancy of the edge state is discontinuous. Furthermore we present numerical calculations on 
a model Hamiltonian to provide additional support for our analytic arguments. 
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In 1988 Haldane pointed out that an insulating crystal 
with broken time-reversal symmetry may exhibit a quan- 
tized Hall conductance even in the absence of a macro- 
scopic magnetic field [1]. We shall refer to such a ma- 
terial as a "Chern insulator" (CI) because it necessarily 
would have a non-zero Chern invariant associated with 
its manifold of occupied Bloch states [2, 3]. While no 
CI has yet been discovered experimentally, there appears 
to be no reason why one could not exist, and theoretical 
models that behave as CIs are not difficult to construct. 
It seems plausible that the current blossoming of interest 
in exotic non-coUinear magnets and multiferroics could 
yield an experimental example before long. 

CIs occupy a middle ground between metals and or- 
dinary insulators. Like metals, their conductivity tensor 
aaf} is non-zero, their surfaces are metallic (as a result of 
topological edge states crossing the Fermi energy), and 
it is impossible to construct exponentially localized Wan- 
nicr functions (WFs) for them [4]. On the other hand, 
only the off-diagonal (dissipationless) elements of aap can 
be non-zero, the chiral edge states decay exponentially 
into the bulk, the one-particle density matrix decays ex- 
ponentially in the interior [5], and the localization mea- 
sure ill [6, 7] is finite [5] as in other insulators. Overall it 
appears natural to regard a CI as an unusual species of 
insulator, but many aspects of its behavior remain open 
to investigation. 

As is well known, the electric polarization P is not well- 
defined in a metal. For an ordinary insulator, its defini- 
tion alternatively in terms of Berry phases or WFs is by 
now well established [8, 9, 10]. For a CI, the absence of a 
Wannier representation removes the possibility of using it 
to define the polarization, and we shall show below that 
there is a fundamental difficulty with the Berry-phase 
definition as well. In view of the presence of dissipa- 
tionless currents and metallic edge states, one might be 
tempted to conclude that P is not well-defined at all in a 
CI. On the other hand, ili is related to the fluctuations 
of P [11], and the finiteness of this quantity [5] suggests 
that the polarization might be well-defined after all. 

The purpose of this Letter is to discuss whether, and 



in what sense, a definition of electric polarization is pos- 
sible in a CI. We demonstrate that the usual Berry-phase 
definition does remain viable if it is interpreted with care 
when connecting it to observables such as the internal 
current that flows in response to an adiabatic change of 
the crystal Hamiltonian, or to the surface charge at the 
edge of a bounded sample. 

For the remainder of this Letter we restrict ourselves 
to the case of a 2D crystalline insulator having a single 
isolated occupied band. The generalization to the case 
of a 3D multiband insulator is not difficult, but would 
complicate the presentation. We also restrict ourselves to 
a single-particle Hamiltonian, noting that the principal 
difficulties in understanding CIs occur already at the one- 
particle level. The lattice vectors ai and a2 are related 
to the reciprocal lattice vectors bi and b2 in the usual 
way (bi • SLj = 2TrSij) and the cell area is S' = |ai x a2|. 

The Berry-phase expression for the electric polariza- 
tion can be written as 

P[ko]-7AlIm/ rfk(uk|Vk|?/k) (1) 

where e is the charge quantum (e > 0), |uk) are the cell- 
periodic Bloch functions, and [ko] indicates the parallelo- 
gram reciprocal-space unit cell with origin at ko (that is, 
with vertices ko, ko-l-bi, ko-l-bi -|-b2, and ko-|-b2). In an 
ordinary insulator one insists on a smooth and periodic 
choice of gauge (relative phases of the juk)) in Eq. (1), 
and P is well defined (modulo eR/S, where R is a lattice 
vector [8]) independent of ko. However, in a CI such a 
gauge choice is no longer possible. To see this, we de- 
compose P[ko] = -Piai -I- P2a2, k = fcibi + /c2b2, and 
ko = Kihi + K2^2, and rewrite Eq. (1) as 

Olih) = -Im dki{uki,k2\dki\uki,k2) ■ (3) 

Eq. (3) is a Berry phase and is gauge independent mod- 
ulo 2tt (independent of ni). This allows us to make an 
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arbitrary choice of branch for di{k2 = K2) and to insist, 
as part of the definition of p['^^\ that 9i{k2) should re- 
main continuous as fc2 is increased from K2 to K2 + 1. 
Since states at (fci, K2) a-nd (fci, K2 + 1) are equivalent, it 
follows that 
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where C is an integer. In fact C just defines the Chern 
number, and the insulator is a CI if C 7^ 0. For simplicity 
we focus henceforth on a CI having C = ±1. 

Using Eqs. (2-3) and similar equations for P2, we have 
arrived at a definition P[ko] that is well-defined, modulo 
eR/S* as usual, even for a CI. However, as illustrated in 
Fig. 1(a), 
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where z is the unit vector along ai x a2 . This dependence 
on ko clearly presents a problem for the interpretation of 
Eq. (2) as a "physical" polarization in the case of a CI. 

However, let us recall how the concept of polarization 
is used. For a normal insulator at least [8], the change of 
polarization during an adiabatic change of some internal 
parameter of the system from time ti to tf is given by 



diJ(t)=Pj)] (modulo eR/5) , (6) 
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where J(t) is the cell-averaged adiabatic current flowing 
in the bulk. A related statement, connected with the re- 
quirement that the charge pumped to the surface must 
be consistent with Eq. (6), is that the charge on an insu- 
lating surface normal to reciprocal vector bi is [9] 



a ~V hi (modulo 6/02) . 



(7) 



Eqs. (6) and (7) embody the attributes of a useful defi- 
nition of P. In the remainder of this Letter, we demon- 
strate that a generalized definition of P, having simi- 
lar attributes, can be given in a CI. We first show that 
Eq. (6) remains correct, provided that the same ko (i.e., 
the same reciprocal-space cell) is used for P*^') and P*-^^ 
in Eq. (6). We also show that Eq. (7) must be modified 
and explain how. We provide numerical tests as well as 
analytic arguments for both claims. 

We begin by giving two arguments for the correctness 
of Eq. (6) in the CI case. First, it is straightforward to 
see that the the contribution to Ji(t) can be computed 
independently for each k2 [9], with the problem in (fci,t) 
space effectively corresponding to that of an ordinary ID 
crystal. Thus, the derivation of Eq. (6) given in Ref. [8] 
goes through unchanged for the CI case. Second, we 
note that the expected result is obtained for the special 
case that the parameter of interest is a spatially uniform 
but time-dependent vector potential A(t). Since a slow 



(a) 

27r 



1 


1 




V 


A/c 

















K2 



K2 + 1 




FIG. 1: (Color online) (a) Sketch of 61(^2) in a Chern insu- 
lator (C = -1-1)- Solid black and dashed red frames indicate 
reciprocal-cell origin chosen at K2 and K2 -|- Afc respectively. 
Dotted lines indicate corresponding averages, proportional to 
Pi . (b) Computed Pi (a) and Vi (a) for the modified Haldane 
model, in units of —0.016/5, for adiabatic (dashed lines) and 
thermal (solid line and symbols) filling. See text. 



turning on of A(t) causes state wue*'' ' 
u\^+(e/hc)Ae'^''' , it follows that 



to evolve into 
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But a time varying vector potential generates an electric 
field E = {-l/c)dA/dt, so that J = {Ce^/h)z x E. The 
transverse conductivity axy is thus quantized in units of 
e^//i, expressing the fact that a CI is a realization of the 
integer quantum Hall effect [1]. 

We further confirm the validity of Eq. (6) by numer- 
ically testing our prediction on the Haldane model [1], 
a tight-binding model for spinless electrons on a honey- 
comb lattice at half filling with staggered site energies 
and complex second-neighbor hoppings chosen so that 
C~l. Using the notation of Ref. [1], we adopt param- 
eters h = 1, i2 = 1/3, = 7r/4, A = 2/3 and the 
lattice vectors ai = aQ{^/3x + y)/2 and a2 = aoy (so 
that oi = 02 = io)- Furthermore, we modify the first- 
neighbor hopping ti —fti{l + a) on the bonds parallel to 
ai -|- a2 so as to break the threefold rotational symmetry 
and allow an adiabatic current to flow as a is varied. The 
compensating ionic charge is assumed to sit on the site 
with lower site energy. 

We consider an infinite strip of the Haldane model A^i 
cells wide and extending to ±00 along y, as sketched in 
the inset of Fig. 2. States ^^^^(r) are labeled by fc2, 
which remains a good quantum number, and an addi- 
tional index n = l,...,2Ni. The dipole moment across 
the strip, per unit length, is 
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where position vector r is decomposed as r = riai -|-r2a2 
and N{k2) is the set of occupied states to be discussed 
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shortly. In the limit of large Ni, we associate the in- 
tegrated current that flows along x in the interior of the 
strip during an adiabatic evolution from a = ai to a = af 
with the corresponding change in Vi , since by continuity 
the charge must arrive at the surface. We then compare 
this with the change of Pi evaluated using a single bulk 
unit cell via Eqs. (2-3) to validate the theory. 

There is a subtlety, however. Neutrality implies that 
N(k2) contains Ni states, but which ones? The problem 
arises because a CI is topologically required to have chiral 
metallic edge states. Our ribbon of CI therefore has one 
band of edge states along its left (L) edge and one along 
its right (R) edge (see inset of Fig. 2). For any given 
a, let fcj {a) be the value of k2 at which L-edge and R- 
edge bands cross. A thermalized filling of the edge states 
would correspond to the thick black curve for case Ofi in 
Fig. 2, where the A^i lowest-energy states are occupied 
at each k2 and ep ^(^2 )• Defining fcj to be the point 
at which the occupation switches between L and R edge 
states, we have fcj = A:^ for the thermalized case. 

In general fcj ip) varies with a. However, fcj cannot 
change during an adiabatic evolution. Because we want 
to "measure" the polarization by the charge that accu- 
mulates at the surface, we specify that the adiabatic evo- 
lution is fast compared to the tunneling time between 
edge states but slow compared to all other processes, so 
that electrons cannot scatter between edges. Thus if we 
thermalize the system at ai and then adiabatically carry 
the system from to af , we arrive at the adiabatic filling 
illustrated by the thick red curve for case at in Fig. 2. 

We thus expect that the change in polarization calcu- 
lated from the right-hand side of Eq. (6) from the bulk 
bandstructure using Eqs. (1-3) should match that given 
by the change of Eq. (9) only if the adiabatic filling is 
maintained. We have confirmed this numerically for our 
modified Haldane model. The polarization as a function 
of a calculated using Eq. (9) and using the right-hand 
side of Eq. (6) is indicated in Fig. 1(b) with black and 
blue dashed lines respectively [12]. Eqs. (2-3) were evalu- 
ated on a 300 x 300 k-point mesh. Eq. (9) was calculated 
using five values of iVi e [25, 70] and then extrapolating 
to infinity, while the k2 integral was discretized with 5000 
k-points. While there is a vertical offset between these 
curves that depends on the choice of kg in Eq. (6), the 
differences APi between different a are correct at the 
level of 10~^. On the other hand, the results obtained 
with the thermalized filling in Eq. (9), shown by the solid 
line in Fig. 1(b), are drastically different. These results 
confirm that the appropriate comparison is with the adi- 
abatic filling, and provide numerical confirmation that 
Eq. (6) is indeed satisfied even in a CI. 

We now turn to Eq. (7). A naive generalization to the 
CI case might be that cr = P[ko] ' bi (modulo 6/02), but 
this cannot be correct. First, the left-hand side should be 
independent of ko, but the right-hand side is not. Second, 
the usual proof for ordinary insulators of the connection 




FIG. 2: (Color online) Sketch of a band structure of a finite 
ribbon of a Chern insulator. Solid regions indicate projected 
bulk bands; thin solid lines are edge states. Black and red 
correspond to a = and a = at respectively; corresponding 
values of fc^ are indicated. Thick lines indicate filling of edge 
states as dictated by fcj , chosen to illustrate system thermal- 
ized at Qi and then carried adiabatically to Off. Inset: edge 
states associated with left (green) and right (blue) surfaces. 



between surface charge and bulk polarization assumes 
that the surface is insulating, with the Fermi level ly- 
ing in a gap common to both the bulk and surface [9]. 
When chiral edge states arc present, the surfaces cannot 
be insulating, so the usual conditions arc violated. 

To show how Eq. (7) can be corrected for the case of 
a CI, let us again consider our Haldane-model ribbon 
at some fixed a. Its surface charge cr can be calculated 
from a = T ■ hi = {S/a2)Vi with Vi evaluated using 
Eq. (9), but its value will depend on the the choice of 
the ^2 at which the occupation of the edge state has its 
discontinuity, so that 

— e 

— dfc2 (V^ifcakllV'nfcs) , (10) 

Jo „^ 



N- 



where J\f is the set of A^i occupied states at /c2 given the 
specified A;| (i.e., the choice whether the L or R edge state 
is included in Af flips as ^2 passes through fcj). 

Since the surface charge theorem of Eq. (7) for ordi- 
nary insulators was demonstrated via the Wannier repre- 
sentation [9], we take the same approach here. However, 
well- localized bulk WFs do not exist in a CI [4] , so we fo- 
cus instead on "hybrid Wannier functions" (HWFs) [13] 
in which the Fourier transform from Bloch functions is 
carried out in the ri direction only. Thus k2 remains a 
good quantum number and the HWF 

WkA-rur2)^^/N'i f dk^ 'FkMri,r2) (11) 

is well localized only in the ai direction. Using these we 
can represent the polarization 
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in terms of the HWF center pjf/' = {Wk^ri\Wk^). We 
require p to be a continuous function of /c2 G [k2, K2 + 1] 
so as to guarantee a result that is equivalent to Eqs. (1-3). 

To make the connection between Eqs. (10) and (12), we 
recast the former by constructing Wannier-like functions 
along the ai direction for the finite-width strip, starting 
from the A^i X A^i matrix 7?,|^^J ,^^ = (V'mfcakilV'nfcs), where 
m,rt g A/" as specified by kl^. The A^i eigenvectors of 

'R}^^^ correspond to states that are Bloch-like along r2 
but localized along ri , which we refer to as ribbon HWFs, 
and the eigenvalues g^^^J locate their centers of charge. 
Using the basis-independence of the trace, Eq. (10) can 
now be rewritten as 
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The similarity between Eqs. (12) and (13) suggests 
that these can be connected. Since k2 is a good quan- 
tum number, each k2 can be treated independently. For 
each k2 we can compare the infinite (bulk) ID system 
described by Eq. (12) with the finite (ribbon) ID system 
described by Eq. (13). The essential observation is that, 
in the limit of large A^i, the HWF centers gjk2 deep inside 
the ribbon converge to the bulk pfca, modulo an integer 
[9]. This is illustrated in Fig. (3), where both sets of 
HWF centers are plotted as a function of fc2 for a ribbon 
of width Ni = 6. Furthermore, the fact that the occupa- 
tion of edge states switches between L and R edge at k2 
is reflected in the discontinuity of ribbon HWF centers 
Qjk^ at fcj. On the other hand, the bulk HWF centers 
are chosen to be continuous across k^- We can account 
for this discrepancy either by including a correction term 
proportional to (fc^ — K2), 
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SpI"'' -heC(fc^ -K2) (mode/a2), (14) 



or by realizing that by the virtue of Eq. (5) this is equiv- 
alent to shifting the reciprocal space origin to 



0-2 



(15) 



as can be seen from the dashed frame in Fig. 3. Eq. (14) 
or (15) is the appropriate generahzation of the surface 
charge theorem, Eq. (7), to the case of a CI, and should 
be correct in large A^i limit for both thermalized and adi- 
abatic fillings as long as the appropriate k2 is used. 

We have also tested the correctness of this formula us- 
ing our numerical calculations on the modified Haldane 
model. Recall that the solid curve in Fig. 1(b) repre- 
sents the surface charge as computed from Eq. (9) for 
the thermalized case. For each a, we first locate fcj us- 
ing 1000 k-points on a ribbon of width Ni = 70 and 
evaluate Eq. (15) with k2 = fc2 using Eqs. (2-3) on a 
250 X 250 k-point mesh. The resulting values are plotted 
as blue dots in Fig. 1(b). The agreement is excellent. 
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FIG. 3: (Color online) Black dots show ribbon HWF centers 

^fk2 ^^"^ ^^'-^ triangles bulk HWF centers pj!^^' and its periodic 
images as a function of k2. Dashed frame corresponds to 

[fe*- 



(^23) choice of origin at discontinuity in g^,^^ , fcj 



In summary, we have generalized the Berry-phase con- 
cept of polarization to the case of a Chern insulator. 
The integrated current flow during adiabatic evolution 
is given by Eq. (6), where the reciprocal-space cell must 
be the same in both terms on the right-hand side. The 
surface charge at an edge of a bounded sample is given by 
Eq. (15), where ^2 specifies the wavevector at which the 
occupation discontinuity occurs in the chiral edge state. 
These results may be of use in understanding the physical 
properties of these topological insulators, and perhaps in 
searching for experimental realizations. 
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